Boundary one-point function, Casimir energy and boundary 
state formalism in D + 1 dimensional QFT 



Z. Bajnok^, L. Palla^, and G. Takacs^ 
February 2, 2008 



^ Theoretical Physics Research Group, Hungarian Academy of Sciences, 
1117 Budapest, Pdzmdny Peter setdny 1/A, Hungary 
Institute for Theoretical Physics, Eotvos University, 
1117 Budapest, Pdzmdny Peter setdny 1/A, Hungary 



Abstract 

We consider quantum field theories with boundary on a codimension one hyperplane. 
Using 1 + 1 dimensional examples, we clarify the relation between three parameters char- 
acterising one-point functions, finite size corrections to the ground state energy and the 
singularity structure of scattering amplitudes, respectively. We then develop the formal- 
ism of boundary states in general D + I spacetime dimensions and relate the cluster 
expansion of the boundary state to the correlation functions using reduction formulae. 
This allows us to derive the cluster expansion in terms of the boundary scattering ampli- 
tudes, and to give a derivation of the conjectured relations between the parameters in 1 -I- 1 
dimensions, and their generalization to Z? -I- 1 dimensions. We use these results to express 
the large volume asymptotics of the Casimir effect in terms of the one-point functions 
or alternatively the singularity structure of the one-particle refiection factor, and for the 
case of vanishing one-particle couplings we give a complete proof of our previous result 
for the leading behaviour. 



1 Introduction 

In this paper we treat quantum field theories with a boundary, which is supposed to be a 
hyperplane (of codimension 1) in flat spacetime. We also treat finite size effects when two 
such boundaries are parallel to each other and separated by a distance L (a.k.a. the Casimir 
effect in a planar situation). 

Such theories have been extensively studied in 1 -|- 1 dimensions since the seminal paper of 
Ghoshal and Zamolodchikov [T], which introduced some fundamental ideas into the subject. 
They defined and characterized the so-called boundary state in integrable field theorie£], which 
makes possible the description of boundary phenomena using the Hamiltonian formalism of 
the bulk theory. Using this concept, they set up the necessary notions to study integrable 

^In Section 3 we give a definition for the boundary state in any boundary quantum field theory in arbitrary 
number of spacetime dimensions, but the details are not necessary for the time being. 
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boundary QFTs by extending the ideas of analytic S matrix theory, factorization and bootstrap 
already well-known in the bulk situation. 

Boundary field theories in 1 + 1 dimensions are relevant to numerous condensed matter 
phenomena (the most prominent examples are the so-called "impurity" problems), to nonper- 
turbative aspects of string theory ("branes") and also are of interest from the theoretical point 
of view as a tool to understand quantum field theory from a different perspective. However, 
boundary problems are also relevant in higher dimensional quantum field theories as descrip- 
tions of surface critical phenomena, and, as mentioned above, in the context of the Casimir 
effect. 

One of the very interesting phenomena that can occur in the presence of a boundary is 
the existence of nontrivial ground state configurations, which approach the bulk vacuum only 
asymptotically, far away from the boundary. In a theory with a mass gap m in the bullo, 
any vacuum expectation value (VEV) of a field cp is bound to approach the bulk value 0o 
exponentially 

(0(t,x))a~</'O + 5ae-"^^ (1.1) 

as a function of the distance d from the boundary, where Qa is a quantity characterising the 
boundary condition labelled by a. We mention that g itself is a physically relevant quantity in 
certain applications, e.g. in the context of thermal Coulomb plasma in contact with an ideal 
conductor electrode where it corresponds to the screened (or 'renormalized') surface charge 
density [2]. The existence of the g term in (II. Ih can be related to the presence of a one-particle 
term in the boundary state of the form [3j 

l^a) = (l + 5a^+(0) + ...)|0) (1.2) 

where A+ (k) - with the normalization [A{k), A+{k')] = (27r)^ y^k^ + m'^6{k - k') - denotes 
the operator that creates an asymptotic particle of momentum k, and (ja is just another 
quantity characterising the strength of the one-particle contribution to the boundary state. It 
is straightforward to derive (see subsection 3.2) the relation 

9a = (1.3) 

where Z is the wavefunction renormalization constant of (j), considered as an interpolating 
field for the bulk asymptotic multi-particle states: 

(O|0(O)|A(k = O)) = y| 

where |0) is the bulk vacuum state and |^(k = 0)) is the asymptotic one-particle state con- 
taining a single particle with the lowest mass m and momentum zero. Note that Z is a bulk 
quantity which is independent of the boundary condition. 

It was also shown in [4j that for 1 -|- 1 dimensional QFTs the presence of a one-particle 
term in the boundary state leads to the following asymptotics of the ground state (Casimir) 
energy on a strip with boundary conditions a and /3 

E^p{L) ~ -mg^gpe-"''' (1.4) 



^In the following we always suppose that this is the case, i.e. we consider only theories with a non-vanishing 
mass gap, for which the formalism of asymptotic states is well defined. 
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while in the absence of such one-particle coupling the energy is expected to decay as e~ 
(which was already known for the integrable case from studies of the thermodynamical Bethe 
Ansatz [5j). In subsection 3.5 we show that this is also valid in general D + 1 spacetime 
dimensions, and provides a more convenient way of defining the quantity ga than the expansion 
(11.21) . The definition of Qa from the one-particle term in the boundary state depends on the 
convention of normalising the one-particle state, while the Casimir energy (|1.4p is a directly 
measurable quantity independent of any conventions in the field theoretic formalism. 

It was already noted in [1] that there is another manifestation of a one-particle coupling to 
the boundary, namely that the one-particle reflection factor has a pole at a special location. 
Introducing the rapidity parametrization of the 1 -|- 1 dimensional energy-momentum (e, p) as 
usual 



it can be argued that the one-particle reflection factor off the boundary (which is the amplitude 
for the process involving a single particle both in the initial and in the final state) must have 
a pole at 9 = with the residue denoted by 



which defines another characteristic quantity ga (which is real for unitary theoriesjj. Ghoshal 
and Zamolodchikov identified ga with ga, but later Dorey et al. [3], investigating one-point 
functions in the so-called scaling Lee- Yang model found numerically the relation 



In the previous papers [U [7l [E] we presented evidence (both analytic and numerical) that this 
relation extends to all 1 -|- 1 dimensional integrable quantum field theories. Strictly speaking, 
however, there existed no field theoretic derivation or proof of this relation up to now. It was 
not clear either whether this relation can be extended to general boundary QFTs. 

The knowledge of the relation between the g parameters has another interesting appli- 
cation. Namely, the classical vacuum expectation value and so the classical limit of ga can 
be calculated explicitely by solving the classical field equations. Relating this quantity to 
the residue of the reflection factor (ga) makes it possible to identify the boundary condition 
corresponding to reflection factors found by solving the boundary bootstrap conditions, which 
is the fundamental idea underlying the work by Fateev and Onofri in [9]. 

The central aim of the present work is to extend these results already known in 1 -|- 1 
dimensional integrable QFTs, to nonintegrable theories and further to boundary QFTs defined 
in any spacetime dimensions. Even in the case of integrable theories, the derivation of (|1.6p 
presented here is the first one that is truly general and starts from "first principles" (i.e. does 
not depend on additional assumptions or approximations). 

We start in section 2 by reviewing the status in 1 -|- 1 dimensional QFTs. We present 
arguments in the integrable case, further to those already made in [1], and then extend our 
considerations to the nonintegrable case using semiclassical techniques, where we show that 
both (|1.3p and (II. 6p can be extended to this case. However, there is no way to tackle general 
nonintegrable quantum field theories without further theoretical developments. 

''Note that the above g quantities have nothing to do with the boundary entropy introduced in [6] which is 
also usually denoted by g and is often called the g-function. 



e = m cosh 6 



p = m sinh 6 



Res Ra (^) = i— 



(1.5) 




(1.6) 
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In section 3, therefore, we develop the necessary tools to address the problem, using only 
general concepts of quantum field theory which are valid in any spacetime dimensions. In 
order to develop the boundary state formalism, in the appendix we derive a set of reduction 
formulae relating the matrix elements of the boundary state with bulk asymptotic many- 
particle states to correlation functions. Using them we can describe the boundary state in 
terms of (analytic continuation of) the boundary scattering amplitudes, analogously to what 
Ghoshal and Zamolodchikov did in the case of 1 + 1 dimensional integrable field theories 
[1]. We present the detailed derivation of the one-particle term in the boundary state which 
leads us to the relation (|1.3p . and of the two-particle term, which is given in terms of the 
one-particle reflection factor. We then use a cluster argument to relate the singularity of the 
reflection factor to the one-point function, i.e. g to g. For the case of 1 -|- 1 dimensions this 
gives a proof of relation (II. 6p "from first principles", valid for integrable and nonintegrable 
QFTs as well, while in the case of D -|- 1 dimensional theories it allows us to characterize the 
nature and strength of the appropriate singularity of the reflection factor. In subsection 3.5 we 
give the derivation underlying the claim in our previous paper [10| where we gave a universal 
expression for the leading asymptotics of the Casimir energy (in the planar situation) valid 
also for nontrivially interacting theories, and also show how these results are affected by the 
presence of one-particle couplings to the boundary. We draw our final conclusions in section 
4. 



2 Examples 

In this section we consider several examples as well as some heuristic derivations (under various 
assumptions) of the conjectured relations among the various g-s. Some of the examples and 
derivations are worked out in specific models, while some of the considerations are model 
independent. Part of the material in this section is on the level of (semi) classical considerations 
while the other is on the level of full QFT. All the examples are worked out in 1-1-1 dimensions; 
they intended to provide a clear understanding of the underlying physics, and motivation for 
the derivations presented in the section 3, which are model independent and valid for generic 
number of spacetime dimensions. 

2.1 Quantum theory examples 

2.1.1 Boundary sine-Gordon model with Dirichlet boundary conditions 

Our first example is the boundary sine-Gordon model with Dirichlet boundary conditions. In 
this model the bulk Lagrangian is written as 

CsG{x,t) = ^ {d^^f - ^ (1 - ^°^^^) (2.1) 

where /3 is the bulk coupling, m is the mass parameter and the field satisfies the boundary 
conditions 

$(X, i)U=0 = K > i)\x=L = 

The fundamental range of the parameters ^ is given by < = fo,L ^ f ■ the 

regime < Att the lightest particle in the spectrum is the first breather, and it has the 
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following reflection amplitude off the boundary [TT] 



p(i).p^ (^) (2^ + 1) i^x-h) sinh(| + z^) 



Here is the rapidity of -B^ while 

Stt 



A = ^ - 1 , and r?o,L = =f(A + l)¥?o,L 



are the parameters characterising the bootstrap solution [T]. Note that R^^\6) has a first 
order pole at = i| originating from the (^) factor. This results in the following coupling 
strength between the first breather and the boundary: 



/ 1 + cos — sin r? , , 

9iiv) = ^\ IT irtan-^ 2.3 

^^'^ Y 1 - cos ^ + sin ^ 2A ^ ' 

(The expression under the square root is always positive as long as A > 1 which is necessary 
for the first breather to exist in the spectrum). However, some care must be taken, because 
the sign of the coupling gi must be opposite on the two ends of the strip since these are related 
by a spatial reflection under which the first breather is odd. As a result, eqn. (II. 4p and eqn. 
(|1.6|) predict that the leading finite size correction to the ground state energy on the strip is 
given by 

Ei^o,^,)=m,ilM^e-' + ... (2.4) 

In [4] this expression was checked by comparing it to the exact ground state energy at large but 
finite L-s. The exact ground state energy of the boundary sine-Gordon model with Dirichlet 
boundary conditions was computed numerically (to high precision) from the NLIE proposed 
recently [l2] and the detailed numerical comparison showed an excellent agreement. This 
comparison was later extended to more general boundary conditions in [13| . 



2.1.2 BTBA in the infrared limit 

In this section we calculate the infrared limit of boundary thermodynamic Bethe Ansatz 
(BTBA) to support the relation (|1.6|) . Although the original derivation of BTBA given by 
Leclair et al. [5] is only valid for the case when there is no one-particle term in the boundary 
state, it can be easily argued that the presence of one-particle terms makes no difference to the 
end result, and this is also supported by numerical studies using comparison with truncated 
conformal space (TCS) by Dorey et al. in pH]. 

In a theory with a single particle of mass m on a strip of length L the BTBA equation for 
the pseudo energy e{9) takes the form 

/oo jn/ , s 

— $ {9 - 6') log (1 + X {0') e~'(^') ) (2.5) 

where ^{9) = —i-^logS {9) is the derivative of the phase of the two-particle ^-matrix S{9)^ 
and X (0) = Ka {9) Kp {9) where 

Ko.{9) = Ro. (4"^) 
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and Ka{0) = Ka{—9). Once e{9) is obtained, the energy of the ground state is given by 

/CO jn , 
— cosh log (l + x(^)e-^('') 

In a theory with nonzero one-particle coupling on both boundaries x has a second-order pole 
at = 0. The logarithmic terms in (|2.5|) and in the ground state energy remain integrable 
and the BTBA makes perfectly good sense, but to obtain the correct asymptotic {L — > oo) 
expression one needs to be careful. For large L we finally obtained in [4] 



oo 



Eaf3[L) = -m——e -m / — cosh ^ log (2.6) 

The remaining integral is an expression of order e"^*"^, while the leading term agrees with 
(|1.4|) (for = gp/2) if ga5/3 > 0. It was already noted in [H] (using a comparison with 
truncated conformal space method) that the BTBA equation only gives the correct ground 
state energy in this case. 

Here we extend this result to values of boundary parameters such that gadfi < 0, applying 
a suitable analytic continuation. Note that in the definition of the boundary coupling gp a 
branch of the square root function must be chosen. In all known cases (e.g. Lee- Yang in 
[T4] or sine-Gordon) there exists a branch choice such that the boundary couplings depend 
analytically on the boundary parameters. In this case a straightforward analytic continuation 
of the TBA equations (following the ideas originally proposed in [15]) from the domain where 
daOp > to the one with gaOp < gives the ground state energy as 

— cosh 9 log (^1 + X % e^'*^^^ j 

where now e{9) is the solution of the equation 

e {9) = 2mLcosh9 - log fj^^ + f ' + ^ (^') ^"^''0 ^^'^^ 

and u is determined by the equation 

Kp{iu)Ka{iu)e-'^'''^ = -1 

For L large one can put e(^) = 2m,Lcosh^, yielding sinu ~ u ~ Id^g/sl ^-mL ^ which in the 
energy expression just fiips the sign of the term coming from the integral part. Thus 

i^apyi^j — +... 

in the parameter range where gadp < 0, which is consistent with eqns (|1.6|1.4I) which is the 
result of the cluster expansion ^4]. 



2.1.3 Bethe-Yang equation 

The aim of the following considerations is to emphasize the difference between the poles of 
the refiection amplitude corresponding to boundary bound states and the one appearing at 
9 = i7r/2, and to give an alternative argument for the validity of eqn. p.6p . 
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Suppose that in the half line theory there is a boundary bound state which is described 
by a pole of the reflection amplitude 



t) — IVq 

in the physical strip < iiq < f . The simplest way to determine the leading finite size 
correction to the energy of the bound state is to confine the theory to a strip of width L by 
imposing identical boundary conditions at the ends of the strip and look for purely imaginary 
solutions of the Bethe-Yang equation of the fundamental particl^ 

ii„(6')i?„(0)e^2™^^'"^^^ = 1 (2.8) 

In order to describe the boundary bound states we continue this equation to imaginary values 
of the rapidity as in [17]. Writing 6 = ivQ + i6 and assuming that 5 — > for L ^ oo one 
finds two solutions for 5 such that the corresponding 6-s are in the physical strip, with the 
following large volume asymptotics: 

The energy of these two solutions above the ground state is 

^± ~mcosi;oTmysinwoe-"'^"'^''o (2.9) 

The interpretation of this result is clear: the two energies correspond to the two (approximate) 
wave functions containing the symmetric (resp. antisymmetric) combinations of the identical 
'half line' bound states localized in the vicinity of the two boundaries (cf. |17|). 

For solutions of the Bethe-Yang equation (12. 8|) that for L — > oo go to = the asymptotic 
behaviour is determined by the residue (II. Sp . Furthermore, writing 9 = + i6, only the 
solution with (5 < is inside the physical strip. Therefore in this case we find the following 
result for the energy 

Eby = mcosh^ = m^e-""^ 

In contrast to the case of a boundary bound state discussed above, the other solution is 
nonphysical. Remembering that Eby is the energy of the given state relative to the finite 
volume ground state of the system, it has a different interpretation than in the previous case. 
Since the true ground state is the symmetric combination of the asymptotic ground states 
localized at the left /right boundary, while the excited one is the antisymmetric one [7j, Eby 
can be identified with the difference between the energies of these two states, between which 
one can switch by changing the relative sign between the two ^-s. Using (II. 4p we obtain the 
relation 

AE = 2m~gle-'^^ = Eby = m^e"™^ 

Note that formally this result is just half of the energy difference Ej^ — in (|2.9I) which is 
due to the special status of the i| pole. This is again consistent with eqn. (II. 6|) . 



''This equation neglects vacuum polarization corrections, but one can argue that for large L-s it correctly 
gives the leading contribution - this was checked in [16] for the boundary sine-Gordon model. 
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2.1.4 Connecting the VEV and the boundary state 



Next we consider the connection between the vacuum expectation value ^(0|<I>(a;, t)|0)/3 of a 
field $ which is an interpolating field for the asymptotic particle states (where |0)/3 denotes the 
ground state of the half line x < theory with the boundary condition f3 imposed at x = 0) 
and the coefficient (s) of the expansion of the boundary state, eqn. (11.21) . in the crossed channel. 
Following the ideas presented in [3], this connection is based on the fundamental (defining) 
property of the boundary state, namely on the equivalence of the correlation functions in the 
two channels. Indeed this equivalence applied to the one-point function reads 

0{O\^{x,t)\O)p = {O\^{y,T)\Bp) y = it, x = ir 

Substituting eqn. (|1.2p one finds 

^mix, t)\o)f, = {o\mm + mmo = mpe'^^^ + ■■■ 

where the first term is the v.e.v. of the field in the bulk (which is an x independent constant) 
and the dots stand for the contribution of the multi-particle terms. The matrix element of 
the field in the second term is basically the normalization of the bulk one particle form factor 
in the closed channel 



mmo = ^) = 




(2.10) 



with Z denoting the wavefunction renormalization constant of Exploiting the connection 
between the coordinates in the two channels one finds finally the asymptotic {x — > — oo) 
expression 

^(0|cl>(x,t)|0)^ = (0|cl>(0)|0) + y|5/3e™^ + 0(e2--) (2.II) 
which gives the following relation 




(2.12) 



2.2 (Semi)classical examples 

Most of the examples and considerations presented so far relied on the integrability of the 
underlying model. In order to extend the relations between the various g-s to nonintegrable 
models we now present some further examples, and also develop methods which do not need 
integrability The general framework in this section is that of the (semi) classical approximation 
of the various field theoretical models. 



2.2.1 Connecting the classical VEV and the ground state energy 

We consider a model described by a scalar field on a strip of width L satisfying Dirichlet 
boundary conditions <I'(0,t) = $0 and <I>(L,t) = $i with the bulk (Minkowski) action 

Auik = j dtdx\^{dt^f -]^{d^^f -Um^ (2.13) 

The ground state is a solution of the static classical equation of motion (satisfying also the 
boundary conditions) 

2(5^) -^W = ^' f^W>0 (2.14) 
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where U{^) describes the self interaction of the scalar field and C is (an L dependent) constant 
of integration. For L — > oo this constant vanishes (C 0) and the ground state tends to a 
superposition of two 'half line' solutions determined by the two boundary conditions. This 
ground state has a simple qualitative description if we assume that U (<I>) has a minimum at 
(normalized as U{^^:) = 0) and that ^o,L are in the 'vicinity' of (meaning that no other 
minimum is between and $o or the scalar field starts at $o then increases (decreases) 
towards reaches a maximum (minimum) $i close to then increases/decreases to $z|j. 
It is important to notice that the L dependence of the ground state energy can be determined 
simply from that of the constant C. Indeed, writing 



E{L) 



and exploiting 




+ f/($) 







1 = -CL+ < 


A/ 




*0 'I'L J 



dv^y2{U{v) + C) 



we get 



L 



<i>i #1 

+ 



dv 



^j2{U{v) + C) 



dE 

11 



thus 



E{L) = E, 



C{L)dL 



(2.15) 



Therefore we determine next the large L asymptotics of the integration constant C(L) using 
the qualitative picture of the ground state. If, in the half line theory, far away from the 
boundary, the solutions with the {^l) b.c. have the asymptotic form 

$(x) 506-"^", Hx) SLe-'"" (2.16) 

then for large L-s, in the central part of the strip, which is far away from both boundaries, 
the scalar field is given by (upto subleading exponential corrections) 

$(x) ~ + goe-"^^ + 5ie'"(^-^) 

(Note that the bulk vacuum expectation value appearing in (I2.16P is a parameter of the 
bulk theory; the dependence on the boundary conditions enters only via the coefficients of 
the exponential terms). Expanding the scalar potential in the vicinity of as U{^) = 
^($-^>*)^ + C'((^>-$*)^), and using eqn. (|2.14ll in the central part of the strip, one readily 
finds 

C{L) = -Im^gogLe-"''' 
Using this in eqn. (|2.15|) gives the asymptotic expression of the ground state energy as 



E{L) = Eoo- 2mgogLe 



-rriL 



where E^o is the sum of the energies of the two half line solutions having the asymptotic 
behaviour in (|2.16ll . This gives 

9 

which is consistent with eqn. (|2.12ll since for the normalization of the field given in (|2.13l) the 
classical value of the wave function renormalization constant is Z = 1. 



'This description is valid if $0 and are both smaller or both greater than $* 
then $1 = "I>, and the solution is monotonic. 



if $* is between them, 
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2.2.2 Semiclassical limit of sine-Gordon model 



In this subsection we consider the semiclassical limit of sine-Gordon model on a strip with 
Dirichlet boundary conditions and determine explicitely the three characteristic i^-s in terms 
of the parameters of the model, thus verifying the conjectured relations among them. 

The bulk Lagrangian is written in eqn. (|2.ip and for simplicity we consider first the special 
case when the model is restricted to the negative half line (x < 0) and the field satisfies the 
boundary condition ^{x,t)\x=o = where the parameter is in its fundamental range 
< f = < f ■ In this case the (semi)classical ground sate is given by a 'half soliton 
standing at the location required by the boundary condition 

$bg {x, a+) = I arctan [g^^^^""^)) where e"'""^ = tan ^ 

Recalling the definition of g 

^hg ~ ge"^^ for X ^ -oo 



we obtain from the explicit solution 



5 = |tan^ (2.17) 



To obtain the semiclassical limit of g coming from the reflection amplitude one has to 
determine the appropriate solutions of the differential equation describing the linearized fluc- 
tuations in the standing soliton background. This was done in [7], [l8], and using these wave 
functions one obtains the classical reflection amplitude 

m-iq iq + m cos ifo 

K\q) = , q = msm.h.0 

m + iq iq — m cos (/jq 

Surprisingly, R{q) has a second order pole at = i^, i.e. for ^ ~ it can be written as 

4tan2(^f) 



R{q) 



We can explain this second order pole in the following way. Since in the semiclassical limit 
the elementary field excitations become identical to the first breather we should compare R{q) 
to the (semi)classical limit of the first breather's reflection amplitude (on the ground state |)) 
(12.21) . which has a first order pole at ^ = i^. However, as shown in [T8], in the (semi) classical 
limit A ^ oo the r/ boundary parameter should also be scaled as 

V = ric\{l + A) 

keeping rjc\ finite, which is related to the boundary value of the field via 

In addition R^^\6) also has a pole at ^ = i — ^) coming from the + |) factor in the 
denominator, corresponding to the hulk process of two S^'s fusing into a of zero energy 
(rapidity = i7r/2), which is then absorbed by the boundary. Clearly for /3 ^ ( A ^ oo) 
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6 ^ and one obtains a second order pole. Therefore for sufficiently small (3 (sufficiently 
large A) in the vicinity of ~ the quantum reflection amplitude behaves as 

R\\^ (0) ~ TT^ ^TTTT^ 7^ TT i Ul = fim — = 

1)^^ (6'-if)((9-i(f -ui)) ' A^oo A^oo2A 
In this case the residue of r\}\9) at ^ = if is 



Res dl\9) = -i— 



while in the (semi) classical limit one obtains 



m 



A 



c, = 4^/-2 tan^ ~ 4^/-tan^ (2.18) 



Combining these expressions with the definition of g one finally finds 

/:=Res 4i)(0) = ,^tan2^ 
2 e=q I) ^ ^ /32 2 

i.e. 

,,^^-^~4./^tan4 
2 V TT 2A 

(where we 'removed' the semiclassical limit in writing the last approximate equality) which 
coincides with the semiclassical limit of (|2.3I1 as was found also in [4]. 

In ^[19] the classical ground state energy of the sine-Gordon model satisfying the Dirichlet 
boundary conditions on a strip of width L 

(in the sector of zero topological charge) was found to be 

E{ipo, (fL, L) ^ Eoo ^ tan — tan — e 

in the asymptotic regime I = mL » 1. Here E'oo is the sum of the energies of the two asymp- 
totic static 'half solitons representing the (asymptotic) ground state on the strip. Defining 
the classical counterpart of g in (11.41) as 

E{L)=E^-m~go~gLe-'^^ 

and comparing to the explicit expression gives 

5i = ^tan^, i = 0, L (2.19) 

Combining eqn. (|2.17|2.18p and (|2.19|) one obtains for the semiclassical limit of the various 
g-s: 



11 



This is consistent with eqn. (11.51) and (|1.6p and - since Z — > 1 in the (semi)classical hmit - 
also with eqn. (|2.11l) . 

Similar results were obtained in [9] for the case of general affine Toda field theories (ATFT) 
where the authors use a different normalization for the fields and the particle states than we 
do. To facilitate the comparison, we give the relations for the case a^^^ (sinh-Gordon) here. 
The action in that paper is normalized as 

^ = I d^x I (d.c^f + M (e^"^ + e-v^'"^) } 

which in our notation means that the classical wave function renormalization constant is 
Z = Air. As a result, from (|2.10l) we have 



(O|0(O)|0 = O) = ^2^+0(6^) 
while their one-particle state , |^)fo is normalized according to 

{o\(t>{o)\e = o)fo = + o{b^) 

which means that \6 = 0) = V^\9 = 0)fo- The residue of the reflection factor is parameterized 
in [9] as 



m 



*2 



which means that D{b) can be written in our notation as 

The form of the one-particle term in the boundary state given in [9] 

|i?) = (l + i?(6)^+o(0) + ...)|0) = (l + :^A+(0) + ...)|0) ^ 9 = ^ 

is then consistent with the relation (|1.6fl^ . In addition, from the asymptotics of the classical 
vacuum solution found in [9] the VEV parameter g in the classical limit 6 — > reads 

^classical = T lim^/7^bD{b) = \/27rgclassical 

which is then consistent with (|1.3p . The results of [9] therefore provide a generalization of the 
arguments of the present subsection to the case of any ATFT. 

2.2.3 Semiclassical description of spontaneously broken theory 

In this subsection we consider the spontaneously broken <I>^ theory restricted by Dirichlet 
boundary conditions to a strip of width L. Using semiclassical considerations we express ex- 
plicitely the (semiclassical limit of the) three characteristic g-s in terms of the parameters of the 
model, thus verifying explicitely the conjectured relations among them in this nonintegrable 
case. 



We mention, however, that the two-particle term in the cluster expansion of \B) is not properly normalized 
in [9], but it plays no role in the issue at hand. 
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The bulk Lagrangian of this model can be written as 



dt J 2\dx J 4 \ X 

where A > and m is real. Introducing z = ^ the well known kink (anti-kink) solution 
'standing' at x = a can be written as 

= t&nh[z - zo), zq = —j= 

It is also well known that the mass of the elementary particle in the 'vacuum' sector is /i = 

First consider this model restricted to the half line x < (or x > 0) imposing Dirichlet 
boundary condition at x = 0: 

which makes it possible to determine the explicit form of two of the three g-s. The essential 
observation is that - for — < $o < least - the (semi) classical ground state of 

the model, $bg, is given by a bulk kink/anti-kink standing at the position required by the 
boundary condition 

vn 

^--^tanh(zo) = $o 

This enables us to determine g which is the parameter characterizing the field in the ground 
state: choosing, say, the solution that far away from the boundary at x = tends to the 
bulk ground state, then g is defined as 

depending on whether we restrict the model to the x < or to the x > half line. Using the 
explicit kink/anti-kink solution one obtains readily 

To obtain g which characterizes the singularity of the reflection amplitude, we restrict 
the model to the x < half line and determine the appropriate solutions of the differential 
equation describing the linearized fluctuations in the standing soliton (kink) background. This 
differential equation is the same as in the bulk case; the only remaining task one is to combine 
the left and right moving bulk wave solutions to satisfy the boundary condition at x = 0. In 
the bulk case the solutions in the continuum spectrum (which is relevant here) can be written 

m 

rjq{z) = e*''^(3tanh^(2: — zq) — 1 — — 3iqtanh{z — zq)) 

where g is a real parameter that determines the frequency of the fluctuation through a;^ = 
m?{2 + q'^/2). We look for solutions in the form 
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that satisfy the boundary condition hk{0) = and also for which the coefficient of e*^^ far 
away from the boundary is one. The classical limit of the reflection amplitude is then defined 

by 

lim hk{z) = e'^^ + Ra{k)e~'^^ 

z— » — oo 

and from the explicit form of hk{z) one obtains: 

„ 2 — fc^ — 3iA; Stanh^fz — zo) — 1 — + 3i/ctanh(z — zq) , „ • , /, 

Rc\(k) = -7^ ; o , k = 2smh6' 

"'^ ^ 2 - A;2 + 3ik 3tanh2(z - zq) - 1 - k^ - 3iA;tanh(z - zq) 

with 9 being the usual rapidity parameter (in terms of which the spatial momentum of the 
particle is p = ^sinh^). 

Rc\{k) has a second order pole at = z^, i.e. for ^ ~ it can be written as 



Ra{k) 



(^-if)2 



This second order pole can be explained in the same way as in the sine-Gordon case. We 
assume that in the full quantum reflection amplitude Rq (6) there is a pole close to q in 
such a way, that in the semiclassical limit (i.e. when A ^ 0) it coincides with the first order 
'quantum' pole at 9 = i^: 



A_ 



i^) ~ 77; 7¥T77; TTT ' J™„ ^ 



Indeed in this case the residue of Rq (9) at 9 = is 



Res Rq (9) = -i- (2.23) 



while classically one obtains 



where 



^cl (k) ~ j2 7^ 



A = -12e-^^° (2.24) 



This mechanism can only work if in the bulk theory there is a fusion of appropriate bulk 
particles to generate the extra pole. Fortunately in the bulk spontaneously broken $^ theory 
there is a rather complex semiclassical spectrum of 'approximate breathers' |20j and u can be 
identified with the fusion angle associated to the second lightest particle B2 as a bound state 
of two copies of the lightest particle Bi : 

2mi cos U = 771-2 

Using the semiclassical formulae for mi and 1712 |20j we get 

m2 9 

cos u = > 1 J 

2mi 32 

i.e. 

u=^^ (2.25) 
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Figure 2.1: Qualitative behaviour of the classical ground state solution 
for A — > 0. Substituting (|2.24|2.25|) into the definition of g gives 

= Res i?q (9) = il^e-^^°, i.e. g = -m^e-^^" (2.26) 
2 0=i| A ^A 

To obtain the explicit form of the third g parameter we confine the model to a strip of width 
L by imposing Dirichlet boundary conditions {^{0,t) = and ^{L,t) = $l) and determine 
the classical ground state energy for asymptotically large L-s: 

E{L) = E^-fi~go~gLe~^^ 

(In the A ^ limit the classical energy gives the leading semiclassical contribution). Assuming 
< <I>o , <1*L < the qualitative behaviour of the static ground state solution is depicted in 
the figure [2TT] (a quantitative description can be given using the methods of [7]). 

The ground state is a solution of the static classical equation of motion (I2.14p with U ($) = 
^($2 _ 1^)2 obtain the L dependence of C explicitely we introduce <I>(x) = ■^ip{x), 

^o,i,L = ■^^o,i,L and C = —^d , in terms of which the condition determining d{L) is: 

= F>,„ + F(^(,.)|o, , F(*(„)|„, = / (2.27) 

V 

where F{(j){v)\a) is an elliptic integral of the first kind with modular angle a and amplitude 
(j){v) with 

ll-d . 2 , l + dl-d-v"^ 

sma = \ , sm (p = ; ^ 

Vl + d l-dl + d + v^ 

For L ^ CO one has d ^ Q and F has the asymptotic behaviour 



FW„)|a) = iln5-l„Ji±^ 
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Using this in eqn. (|2.27ll determines the L dependence of d, which finally leads to 



E{L) = Eoo- 8^—— — e 

A l + (fo l + (pL 

The boundary conditions require tanhzo = —ifo and tanhz/, = ip^, leading to 

5o = -m^e2^°, = -m'^e-^'^ (2.28) 

Combining eqn. (I2.22|2.26p and (|2.28ll one obtains again that the semiclassical limits of the 
various g-s satisfy (I2.20p . i.e. 

even in this nonintegrable model. 

In passing we note that one can consider the theory in the symmetric phase (i.e. without 
spontaneous symmetry breaking, obtained formally by changing the sign of in (|2.21|) ) 
along the same lines as we did in the previous subsections for the sine-Gordon and for the 
spontaneously broken theories here. In this case one can determine - using semiclassical 
considerations - the g coming from the VEV of the field in the ground state, as well as 
the g defined from the ground state energy on the strip and verify that they satisfy the 
conjectured relation. For the g determined from the refiection amplitude one encounters the 
following problem: although the differential equation of the small fluctuations around the 
ground state with the right boundary conditions can be solved and the classical reflection 
amplitude obtained, it exhibits the same second order pole at 6* = as we encountered in the 
sine-Gordon and spontaneously broken theories. However, in the symmetric phase of 
theory there is no analogue of the bulk soliton/kink solutions, and thus neither of their breather 
bound states (exact or approximate). Despite this, the 'pole merging' mechanism described 
in the previous cases still exists, but the pole which collides with the one already located at 
6 = in the semiclassical limit comes from outside the physical strip, as appropriate analytic 
continuation of the reflection factors to this regime shows. 



3 Boundary state formalism and proof of g-g-g relations in D + 1 
dimensions 

In this section we analyze quantum field theories in D + 1 dimensions in the presence of a 
D dimensional fiat boundary. The correlators are defined via the path integral approach, 
and are then expressed using two alternative Hamiltonian descriptions of the system. If the 
boundary is in space ('open' channel) the information on the boundary condition is encoded 
in the refiection factors, which corresponds to the on-shell part of the correlators. However, 
when the boundary is situated in time ('closed' channel), the information is contained in the 
boundary state. By deriving a reduction formula in the 'closed' channel we are able to relate 
the matrix elements of this state to the on-shell part of the correlators and thus to the refiection 
factors. The knowledge of the boundary state, especially the one-particle boundary coupling, 
makes it possible to calculate the large distance behaviour of the one-point function, as well 
as to determine the leading finite size corrections of the ground state energy (Casimir effect). 
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Using the clustering property of the two-point function, these two quantities can be related 
to the singular part of the one-particle reflection amplitude. In this section we focus our 
attention on one- and two-point functions; the analysis of the multi-pont correlation functions 
is relegated to the appendix. 



3.1 The concept of the boundary state 

Consider an Euclidean quantum field theory of a scalar field $ defined in & D + 1 dimensional 
half spacetime, parameterized as {x < 0,y,r), in the presence of a codimension one flat 
boundary located at x = 0. The correlation functions defined as 

- N\ f V<^^{xi,yi,fi)...^{xN,yN,rN)e-^^'^^ ,^ _ 

mxi,yi,ri)...^{xN,yN,rN)) = ^ ^ / P^g-^I^l ^ ^ ^ 

contain all information about the theory. The measure in the functional integral is provided 
by the classical action 

Sm = J dr P dy y° dx {^{d.,^f + \{dy^f + \{d^f + U{^)^ + C/b($(x = 0,y,f)) 

which determines also the boundary condition via the boundary potential Ub- 

This Euclidean quantum field theory can be considered as the imaginary time version of 
two different Minkowskian quantum field theories. We can consider y = it as the imaginary 
time and so the boundary is located in space providing nontrivial boundary condition for the 
field We refer to this description as the open channel. The field $ is an operator valued 
distribution that satisfies the equal time commutation relation 

[<l>(x, i, r), , t, r ')] = — X )(5(r — r ') ; x,x <0 

The space of states in this Hamiltonian description is the boundary Hilbert space TIb deter- 
mined by the configurations on the equal time slices. TLb contains multi-particle states and 
is built over the boundary vacuum |0)b by the successive application of the particle creation 
operator^. In the asymptotic past the particles do not interact and behave as free particles 
travelling towards the boundary; thus 

'nB = {4n{kiM)---atn{kN.kNmB , h>--->kN>o} 

where the operator af^{k, k), normalized as 

[ain{k,k),af^{k',k')] = {2TT)^uj{k,k)6{k-k')6{k-k') ; k,k'>0 
creates an asymptotic particle of mass m with transverse momentum k and parallel momen- 



tum k. The corresponding energy is U!{k,k) = \/k'^ + k'^ -\- m? = y k'^ + mcs{k)'^ . In the 
Heisenberg picture the time evolution of the field 



'^One can also introduce particle-like excitations confined to the boundary [21 j . but here we do not consider 
them for simplicity. 
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is generated by the following boundary Hamiltonian 

fO 



Hb = j dr j dx Qn? + + 1(9^)2 + U{^)^ + Ub[Hx = 0)) 



The correlator (13. ip can then be understood as the matrix element 

yi, n) . . . ^{xN, VN, tn)) = B{0\Tt ($(xi, ti, fi) . . . <^>{xn, ^tv, rN)) \0)b 

where Tt denotes time ordering with respect to time t, and the vacuum \0)b is normalized to 
1. 

Alternatively we can consider x = ir as the Minkowskian time. In this case the boundary 
is located in time and we can use the usual infinite volume Hamiltonian description. This is 
referred to as the closed channel. The Hilbert space is the bulk Hilbert space H spanned by 
multi-particle in states 

= {Atni'^uki)...A+{KN,kNm , A:i > ••• > A;jv} 

where the particle creation operators are normalized as 

[Anit^,k), At^{K',k')] = {2TTfu;{K,k)6{K-K)6ik-k') 

Time evolution 

cI>(T,y,f) = e^^^f(0,y,f)e-^- 
is generated by the bulk Hamiltonian 

H = J dfj"^^ dy Qh^ + + + ui<^)^ 

The boundary appears in time as a final state in calculating the correlator (|3.ip : 

($(a;i, yi, n) . . . <^{xn, Vn, tn)) = {B\Tr ($(ri, yi, n) . . . $(t^, y^v, r^)) |0) 

The state {B\ is called the boundary state, which is an element of the bulk Hilbert space and 
is defined by the equality of the two alternative Hamiltonian descriptions 

{B\Tr ($(ti, yi, fi) . . . <^{tn, Vn, vn)) |0) = B(0|rt ti, fi) . . . $(x^, t^, ttv)) |0)b 

where the correspondence is valid if {iT,y) is identified with {x,it). Using asymptotic com- 
pleteness the boundary state can be expanded in the basis of asymptotic in states as 



{B\ = (0| J 1 + K'AUO, 0) + [ ——^^—^R\k, k)A^n{-K, -k)Ain{K, k) + ...\ 

[ Jo J {2Tr)^-^uj{K,k) J 

(3.2) 

which we refer to as the cluster expansion for the boundary state (where due to translational 
invariance only bulk multi-particle states with zero total momentum can appear). The bars on 
top of the K coefficients indicate that the above expansion is that of the conjugate boundary 
state. 
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3.2 One-point function and 

The one-point function of the field, due to unbroken Poincare symmetry along the bound- 
ary, only has a nontrivial dependence on x and can be written (using the 'open' channel 
formulation) as 

Bm{x,t,r)\0)B = Gl,y{x) 
Going over to momentum space by Fourier transformation 

and analyzing the all order perturbative expression the following form is obtained: 

Unbroken translational invariance along the boundary is manifested in the conservation of 
energy and parallel momentum. Summing up the contributions of the outer leg the bulk two- 
point function, a pre-factor G'^^^^ can always be isolated. The remaining factor has a part 
that preserves the transverse momentum (which is the same as the bulk one-point function 
^luik) ^'^^ another one which depends on the boundary condition and violates transverse 
momentum conservation. The bulk two-point function has the well-known Kallen-Lehmann 
representation in terms of the spectral function a{m) 

iZ r°° , > ia(m') 

dm 



Gluiki^, k, k) = — z- + / 

_ _ _ _|_ 



2m U) 



— k"^ — k'^ — m"^ + ie 



where the wave function renormalization constant Z characterizes the strength of the on- 
shell part. Performing the energy-momentum integration and picking up the pole terms the 
one-point function can be expressed as 

Bm{x,t,mB = miOm - ^SL^M^"^^ - l2dm"-^^Bl,^y{im')e^'^ (3.3) 



B{Omx,t,r)\0)B = mmO)+ge"'' ; 5 = -— ^LyM (3-4) 



Prom this we can read off the leading large distance behaviour (dominated by the on-shell 
part): 

iZ 
2m 

Our aim now is to connect the quantity g to the coefficient of the boundary state (|3.2|) . 
This can be accomplished by calculating the matrix element 

{B\Al{k,k)) = {2TTfL0{K,k)6{K)S{k)K' 

in the closed channel. Using the reduction formula derived in the appendix, it can be expressed 
in terms of the one-point function as 

{B\Af^{K,k)) = dy r dr I e-*^('''^)^+*'^^+^^^^ 

-00 J— 00 J 



\dl-dl-B^ + m^ - 5{T){dr + iuj{K, ^))} (S|$(r,y,f)|0) 
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In the appendix it is also shown that only the on-shell part of the correlator contributes to 
the quantities K, and thus one can substitute 

{BMT,y,r)\0) = BiOMx = ir,t = -iy,mB « -^e™^B,Vj/M 
Plugging this expression back we obtain 

= [Gboundaryiim)] = 9 

from which using (|3.4|) we can establish the relation 



9 = y^9 

We note that the derivation remains valid if the Lagrangian field ^ is replaced by any bulk in- 
terpolating field for the asymptotic particles and its appropriate wave function renormalization 
Z. 

3.3 Two-point function and the relation between R[ and K"^ 

The two-point function in the open channel (using the unbroken spacetime symmetries) can 
be written as 



B{0\T{'^{x,t,r)^{x ,t ,r ))\0)b = 

duj f dk f dk f dk' 



iu^{t~t )~ikx~ik X ~ik(r-r )gI.{uJ, k, k',k) 



2tt J (27r)^-i J 2tt J 27r ^'^'^^ 

after integrating out uj' and k' using the delta-functions in (IA.3|) . where 

Gldryi^, k, k',k) = 2Tr5{k - k')Gl^iki^, k, k) + Gl^i^i^, k, k)B^a^y{uj, k, k' ,k)Gl^iki^, k' ,k) 

Its on-shell part determines the two-particle reflection factor, which is the probability ampli- 
tude of a transition from a one-particle initial state \k ,k )in = af^{k, k)\0) b into a one-particle 
final state \k^k)out = ^tutik ,k)\^) b and can be expressed as the matrix element: 

out{k, k\k',k')in = {2TT)^uj{k, k)6{k - k')6{k - k')R\{uj{k, k), k) 

where Rl is the one-particle to one-particle reflection factor (cf. (|A.5j) ). In order to simplify 
notation (and to conform better with the conventions of section 2) we drop the indices and 
denote R\ simply by R from now on. The boundary reduction formula connects it to the 
correlator a^ 

out{k,k\k,k = out{k,k\k,k )out 

-2Z-iy"° dx j dt j dfe-^'^('='^)*+*^^^cos(A;x){a2-52-a^ + m2 + 5(x)a^} 

dx jdtj (ife-^'^('=''^')*'+*^'^^cos(A;'x'){a2, - - + + 6ix' )d^>^ 

B{0\Tt{<!>{x,t,r)^x' ,t ,r'))\0)B 



*The normalization of the creation operators agrees with [10] but differs by a factor \/2 compared to |21) 
and this affects the reduction formulae. 
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As was explained in [21] and is shown in the appendix only the on-shell part contributes. As 
a consequence we can keep the on-shell part from the bulk two-point function and represent 
the correlator equivalently as 



B{0\Ttmx,t,r}^x' ,t ,f'))\0)B ~ 

dujdk dk iZe-^^^^-^V^^^~^'^ 



2tt 2tt (27r)^-i ^2 _ ^2 _ ^2 _ ^2 + 

Plugging this expression into the reduction formula, the bulk part of the two-point function 
cancels the disconnected piece and the reflection factor turns out to be 

-> Z -> 

R{uj{k, k), k) = —Bl^^y{uj{k, k), k, k, k) 

Let us connect this quantity to the quantity in the boundary state (|3.2p . considering the 
following matrix element: 

(k, (k, ')) = (27r)^t^(K, k)5{K + K)5{k + k')K\K, k) 

where k > k is assumed. Applying the result of the reduction formula presented in the 
appendix we have 

(S|A+(N^)a+(^',^')) = 

(2Z)-i P dy j dr j dre~^'^'^'^^^^^+^^y+^''^ \dl - d'^ - P + m'^ - (5(r)(itj(K, k) + 5^)} 

dy' f dT [ e~M'^',fc')-'W|/'+ifc'r' 1^2^ _q2^ _ d'\ _ S{T'){iuj{K ,k') + d^>) 



OO J —CO 



(5|r,($(r,y,f)$(r,y,f )|0) 

Using the on-shell part of the two-point function and exchanging the role of space and time 
as X = ir and y = it we obtain 

(5|r,($(r,y,f)$(T',y',f')|0) = B(0|rt($(x, t, f)$(x', t', f '))|0)b « 



dujdk dk iZe^^^y-yV^^^^^'^ 



27r 27r (27r)^-i ,^2 _ ^2 _ ^2 _ ^2 + 

where additionally we also exchanged the role of the energy and transverse momentum as 
oj ^ —ik. In the reduction formula the bulk part depending on r — r does not contribute 
and we obtain 

- -> Z ^ _ _ 

K'^{K,k) = — ^B^a^y{-iK,iuj{K,k),iuj{K,k),k) 

2iLo{K, k) 

which means that the relation between and R is 

K'^{k, k) = R{uj — > — ZK, k iuj) 
In two spacetime dimensions one can use the rapidity parametrization 

to = m cosh , K = m sinh "i? 

and then 



(^) =Ri^i-+^ 

which is the same as the relation derived by Ghoshal and Zamolodchikov [I]. 
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3.4 The connection between g and the reflection factor 

The singularity property of the two-particle reflection factor generally follows from unitarity, 
similarly to the well-known bulk theory of the analytic S matrix. Unitarity of the reflection 
operator TZ can be expressed in terms of the interaction matrix T as 

7^7^+ = (1 + ir){i - ir+) = i 

Calculating the one-particle matrix element of this relation between states \k, k)B and \k' , k)B 
and inserting a resolution of the identity we have 

2QmT{k,k){2Tr)'^uj{k,k)6{k-k')S{k-k') = B{k,k\r\0)B B{0\T-^\k' ,k')B + 

l2^ i 2^^^B{k,k\T\q,q)BB{q,qlT^\k',k')B + ... 

where the second term on the right hand side corresponds to boundary bound states (degrees of 
freedom propagating along the boundary) which are pole singularities located at A; = iun (j^^ 

with < Un{k) < m. Unitarity (as expressed above) gives the following relation for the residue 
of the reflection factor: 



2Qm 



ojiiu, k)R{iu, 0) = 27r^^^^^^(^(fc - i«)C„ (k^ Cn (fc)^ 



where C„ (^k^ is the strength of the on-shell vertex corresponding to the creation of the 
boundary bound state by the bulk particle. 

Using the unbroken spacetime symmetries we can extract energy and parallel momentum 
conservation which in the first term gives 

B{k,k\T\0)B = B{k,k\n\0)B = {2t^)^ 5(k)5{u{k,k))R^{im) 

This term vanishes at physical energies but we can try to make an analytical continuation 
via the boundary reduction formula and express from the one-point function. The energy 
delta function can be rewritten as 

6iuj{kM = ^^^^^(5(fc - im) 
k 

and the kinematical pre-factor vanishes exactly at the momentum value it is concentrated 
on. This shows that the singularity of the two-particle reflection matrix at A; = im and 
fc = is not a pole in the momentum variable k. We see that in this case the unitarity 
argument does not connect the singular part of the reflection factor to the one-particle emission 
amplitude and so to the on-shell part of the one-point function. The fact that there is an 
essential difference between the behaviour of this particular type of singularity and the ones 
associated to boundary bound states is consistent with their different treatment in the Bethe- 
Yang analysis presented in subsection 2.1.3. Therefore we choose an alternative way to connect 
the strength of the singularity in the reflection factor to the VEV of the fleld, making use of 
the clustering property of the two-point function. 
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In the Euclidean regime the on-shell part of the two-point function, which determines its 
behaviour in the asymptotic regime far from the boundary, reads as 



B{0\^{x,y,r}^x',y'y)\0)j, « (3.5) 



dp dk dk Z 



where we supposed t > t and introduced co = —ip, t = iy. Cluster property implies that for 
large temporal separation the Euclidean two-point function satisfies 

' ^/ ^ . ' 'Jm^s . ^ f -^\y-y \ 



B (0| $(x, y, f)$(x ,y,r) \0) ^ = b (0| $(x, y, r) lO)^ ^ (0| ^'(x , y , r ) lO)^ + O [e 

with some characteristic scale p that corresponds to the gap in the spectrum above the vacuum. 
The presence of the disconnected piece signals the nontrivial vacuum expectation value of the 
field ^> (which is time independent dy (0| ^{x,y,f) |0) = due to y-translational invariance, 
and is also independent of f for a similar reason). In the asymptotic regime we then expect 
the following behaviour 



B (0| $(x, y, fMx , y'y) 10)^ ~ ^^e^^ ^ |y _ y | ^ oo (3.6) 

In a free theory, with Z = 1 and R{k, k) = ±1 (corresponding to Neumann/Dirichlet boundary 
conditions), the integral (|3.5p can be evaluated explicitely with the result 



~2 m(x+x ) 



— {Kq (mr^) ± Kq (mr+)) , r± = \/{y — y'Y + (2: i x'Y + {r — r)"^ 
27r V 

which decays exponentially when |y — y'| — > 00 and so there is no disconnected piece. In fact, 
this remains true as long as R{k, k) is regular at A; = ±.im as will be implied by the analysis 
below. 

Here we have to make an Ansatz for the singularity type of the reflection factor. The 
expectation from its structure is that it provides the needed clustering, moreover it should 
match with the exact results available in two dimensions. In two dimensional integrable 
theories the singularity of the reflection factor at = ^ is a pole in the rapidity variable 
[k = msinh^) of the form 



d — coshi? 



Changing the rapidity variable to the momentum this corresponds to the behaviour 



R{k) 



Observe that the singularity in the variable k is not of a pole type, but is milder, just as 
expected from the unitarity argument. In higher dimensions we expect this singularity to be 
connected with the virtual one-particle emission (just as in the 1-1-1 dimensional case), and 
due to parallel momentum conservation it is expected to have the following form: 

R{k,k) mgy2 (2TT)^5{k) (3.7) 

yk"^ + kP' + m? 
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Such a term introduces a singularity at k'^ = —rr? where the on-shell value of p is which 
means a contribution which is constant in the temporal separation \y — y \. We remark that 
in order to obtain this behaviour from the analytic R-matrix theory the classification and 
analysis of Coleman-Thun diagrams at this particular kinematical point is necessary which is 
outside the scope of the present considerations (for a general exposition see [21|). 

We can calculate the effect of this singularity in the following way. For definiteness, let us 
specify y > y : then the integration contour in (13. 5|) can be closed in the 9m p < half-plane 
giving the result 

Bm{x,y,r)Hx\y\r^m^ = - f ^ ^J -V^iv^v') [_^Jl=,-^K^^^')\ 

J 2v fe^ + L yk^ + J 

where we have also performed the trivial k integration. Remembering now that x + x' < we 
can close the contour in this term in the lower half-plane which gives 

B (0| y, r)^x' , y', r^) |0)^ ~ ^^e™(^+^') + . . . 

o 

and comparing to (13. 6p it follows that 



and so 




It is useful to stress that for D > 1 this result shows just the consistency of the assumption 
(13.71) . In 1 + 1 dimensions, however, the exact reflection factors of many integrable theories 
are known explicitely, and exhibit a pole at = ivr/2, and in this case the cluster argument 
is in fact a field theoretic proof of a relation between g and g (valid also for the nonintegrable 
case), conjectured in [3] and checked using both theoretical and numerical arguments in [4]. 



3.5 Finite size energy correction from the boundary state 

We now calculate the ground-state energy per transverse volume, Eq^{L), of the system 
confined by two hyperplanes to the interval < x < L and subject to boundary conditions 
labelled by a and f3 at the two ends. In doing so we compactify the other directions to circles 
of perimeter R with periodic boundary conditions and calculate the partition function in 
two different ways, corresponding to the two different Hamiltonian descriptions of the system 
introduced previously. In the open channel (where y = —it is the imaginary time) the partition 
function can be written as 

Z{L,R) =rr (6-^^-/3 (^'^)) 

which for large R behaves as 

lim Z{L, R) = e'-^o^C^)^ + small correction 

_R— >oo 

where V = R^ . In the closed channel (where x = it is the Euclidean time) the partition 
function is given by the following matrix element. 

Z(L,i?) = (i?„|e-^^(«)|i?^) 
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Inserting a complete set of eigenstates of the periodic (bulk) Hamiltonian H{R) 

^-^ (n\n) 

We concentrate on the leading finite size correction to the ground state energy Ef{L) and 
take L to be large. As a consequence the low lying energy levels dominate the sum: 



^ ^ {Bo,\K,k){K,k\B,3) ^ 
^ {K,k\K,k) 



where the sum is over one-particle states with momentum (k, fc). The finite volume restricts 



the momentum to be k = %n, and ki = ^rii and the normalization of the creation operators 



becomes 

[Ain{K,k),Af^{K ,k)] = Vu:{K,k)5^^^i5^^^, 
Using the form of the boundary state (13.21) yields 

Z(L, R) = e-^o(^)^ [1 + mVKlKle-'^^ + . . . ] 

We normalize the ground-state energy with periodic boundary condition to zero: Eq{R) = 0. 
As a result the ground state energy per transverse volume at leading order in L has the form 

EfiL) = -mkiKy^^ + ... (3.8) 

If one of the K^-s is zero then the leading correction comes from two-particle states. Although 
we explained this situation in [U [TO] , for completeness we recall the derivation of the leading 
finite size correction. Using the cluster expansion (13.21) . the partition function can be written 
as 



Z{L,R) = e~^°^^^^ 



{Ba\K,k,K ,k){K,k,K ,k |^/3) ^-^(^,fc)L-..(«'.fc')L , 



The spectrum of the possible k, k has to be determined by solving the scattering problem of 
the two particles in volume V exactly. However, in the infinite volume limit the interaction 
between the particles can be neglected. (For large L the main contribution to the ground state 
energy comes from the low lying energy levels, for which this approximation becomes exact 
as i? ^ oo.) Using the explicit form of the boundary state, the result for the Casimir energy 
(per unit transverse area) is 



as we explained in our previous paper [10| . where it was derived using the assumption that the 
modes indexed by k can be treated as independent two-dimensional degrees of freedom (their 
interaction only entering higher order terms explicitly). The present derivation, however, 
completely dispenses with this additional assumption. 
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4 Conclusion 



In this paper we have considered the relation between one-point function, finite size corrections 
to ground state energy and the analytic structure of scattering amplitudes in a general bound- 
ary QFT. Our main result is the derivation "from first principles" of the relations (|1.3|1.6p 
between the parameters characterizing these quantities for a general D + 1 dimensional quan- 
tum field theory. 

Besides giving a solid theoretical background to preexisting results, we developed the 
formalism of the boundary state for the D + 1 dimensional case, together with a new set of 
reduction formulae relating the boundary state to the correlation functions, which can be used 
to express the cluster expansion of the boundary state in terms of the boundary scattering 
amplitudes. 

The resulting formalism can be used to address various issues in boundary quantum field 
theory. In addition to the derivation of the relation between the g parameters and its exten- 
sion to general field theories, we have shown that it can be used to derive the large volume 
asymptotics of the Casimir effect. In particular in the case when there is no one-particle 
coupling to the boundary, we proved that the result of [lO] for the leading behaviour is indeed 
universal, even for interacting fields. We extended this result by giving the leading term in 
the Casimir energy (|3.8p for the case with nonvanishing one-particle coupling, and used our 
results to relate it to the asymptotic behaviour of the vacuum expectation value of the in- 
terpolating field. It is important to note that it relates the planar Casimir effect to physical 
quantities (vacuum expectation values, reflection factor) that are calculable in the infinite vol- 
ume boundary quantum field theory, which is a much simpler setting than the finite volume 
case. In addition, the expression is free of any ultraviolet divergences right from the start, the 
underlying reason being that it involves only relations between physical (i.e. renormalized) 
entities, as was already discussed in [10]. It also noteworthy that it gives a general formula for 
the dependence of the Casimir effect on the material properties of the boundary (characterized 
by the reflection factors) in the planar setting. 

In principle, the cluster expansion for the boundary state (13.2(1 provides a systematic large 
volume expansion of the Casimir energy (and other quantities, such as vacuum expectation 
values) if higher particle terms are included. However, in the case when the reflection factor 
has a singularity corresponding to the one-particle coupling, even the two-particle term is 
divergent in general. This is a sort of infrared divergence, which in the case of the ground 
state energy is known to be eliminated by resummation of the series using e.g. thermodynamic 
Bethe Ansatz (cf. [4j where a regular large volume expansion is derived from boundary TBA), 
but this method only works for integrable theories (or for theories in D + 1 dimensions with 
trivial bulk and completely elastic boundary scattering). It is plausible that something similar 
happens in the case of nonintegrable theories as well, but there is as yet no method to perform 
a resummation of the relevant terms, which remains an interesting open problem. 
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A Reduction formulae 



In this appendix we recall the derivation of the reduction formula in the open channel and 
show its physical meaning, namely that it connects the reflection factors to the on-shell part of 
the boundary Green functions. Then we derive an analogous reduction formula in the closed 
channel, and show that only the on-shell part of the boundary Green function contributes to 
the boundary state, thus a relation between the boundary state and the reflection factor is 
established. 



A.l Reduction formula in the open channel 

The multi-particle reflection matrix is defined as the matrix element 



out 
B 



{pi,Pi; ■ ■ ■ ■,Pk,Pk\atniqi,qi)\q2,q2; ■ ■ -^quqiYB (A.l) 



The creation operator of the initial particle state can be written in terms of the asymptotic 
field as ^ 

o+(A:,fc) = — )= [ dxcos{kx) [ df e*^^^e-^'^('='^)*W^i„(x, t, f) 

\/2 J-oo J 

Using the property that the interpolating (interacting) field approaches the free asymptotic 
fields as 

t, r) 2'^/2$i„/o„i(j;, t, r) for t ^oo 
we obtain the following form of the reflection matrix 

Rk = l3Hpi,Pi,---'^Pk,Pk\atutiqi^qi)\q2,q2-,----,qi,qi)B + 

f I j '^*^t{cos(gix)e-^^(«i'«^i)*WB(ont|$(x,t,f)|m)B} 

where B{out\, (Im)^) is the shorthand form for °b\pi^Pi'^ ■ ■ ■ ;Pk,Pk\ and 1^2,92; • • • ^quqiYSi 
respectively. It is necessary to be careful when performing the partial integration and keep 
the surface term. The connected part turns out to be 

dx / dtdfe-"^^'i^''^^^*e"l^''cos{qix){df -dl-d'^ + m^ + 5{x)d^] B{out\^x, t, f)\in)B 

(A.2) 

Repeating the same procedure one can express the reflection factor as the product of integro- 
differential operators acting on the Green functions. The connected part is related to the 
A^-pont function, while the disconnected one appears when at least one of the incoming mo- 
mentum (say {qi-,qi)) coincides with one of the outgoing momentum (say {pj,Pj)). It contains 
a delta function singularity {27r)^ u>{qi, qi) 6 {qi — Pj)S{qi — Pj) whose coefficient is related to 
the N — 2 point function. Here we would like concentrate on the connected part and to show, 
that the operator in (IA.2j) truncates the leg of the momentum space Green function and puts 
its momentum on-shell. 

The momentum space Green function is deflned by Fourier transformation 

B{0\Ttmxi,h,ri)...^{xN,tN,fN))\0)B = 

i-oo 27r y_ 2vr (2vr)^-i' (^....(i'^o fe., A;J) 
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Inspecting the perturbative expansion we can always write 



disc. + Jl Gl^if,{ujj, kj, kj)B^^^y{{uj^, h^ki}) 



(A.3) 

where the disconnected part contains at least one particle which is unaffected by the boundary, 
that is a term proportional to {27r)^6{ki — kj)6{ki — kj). H Using the Kallen-Lehman form of 
the bulk two-point function the contribution of a single leg can be written as 



dio dk dk 



iujt —ikx —ikf 

e e e 



iZ 



u)'^~k'^ — k'^—mP-+it 



ia-(rn )dm 



2m t^2_fc2_fc2_ 



(A.4) 



B^,,Ju;,k,k,...) 



We are interested in the action of the operator 

r dx j dtdfe-''^^'^''^^e''^cos{qx){d^ -dl-d^ + 6{x)d^} 

appearing in (IA.2P on (IA.4j) . Calculating the derivatives dt, d and performing the i, r integrals 
explicitely, we obtain {27r)^6{io{q, q) — io)S{q — k), and thus the io and k integral replaces every 
io by io(q,q) and every k by q. So it remains to be shown that the operator 

rO 

dxcos{qx) {—9^ — q'^ + Six)dx} 



I 

J —I 



operator truncates the leg 



dk 
— f 
2tt 



-ikx 



iZ 



k^ + ie 



+ 



ia{m )dm 



2m Q 



/t2 + 



BMrv(^{q,Q),k,q,...) 



where q 



'2 



q^ +Tin? 



m 



'2 



The actions of —d'^ — q^ and S{x)dx can be computed separately. 



The second is simpler: it gives a factor —ik and substitutes a; = in the k integral. The term 
— q^ gives a factor k'^ — q^ but now the x integral is nontrivial and we use that the contour 
can be closed on the upper half plane giving 

fO 



dx cos{qx)e 



-ikx+ex 



1 



+ 



-i{k - q) + e -i{k + q) + e 
i [V^ + + TT {6{k -q)+ 5{k + q)) 

\ k—q k+q / 



Consider first the terms containing the spectral density a. Due to the pre-factor fc^ — the 
delta functions do not contribute and the principal value can be replaced by the function itself, 
which just kills the other a term coming from the b{x)dx term. Now we analyze the terms 
containing Z . In the formula for 8{x)dx we write 



—ik 



iZ 



q^ — k"^ + ie 



Z 
2" 



1 



+ 



1 



k + q-\- ie k — q — ie 



^This definition of the disconnected part is slightly different from the one used for the two-point case in our 
previous work [21] but is more convenient for multi-point functions. 
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The terms with the denominator k + q + ie cancel, while the ones with denominators k — q±ie 
combine together to give delta functions resulting in 

which, combined with the pre- factor '^^J^ in the reduction formula, gives the contribution of 
one leg as 

Reduction of an outgoing particle gives rise to the same effect, so collecting the contributions 
of all legs the connected part is simply 

i j i j 

BMry{^{Pi^Pi)^Vi,Pi,^{<lj,qj),(lj,%) (A. 5) 

In addition to the connected part treated above, there are two types of disconnected parts: one 
from the reduction formula itself, and the other originating from the disconnected part of the 
Green function (|A.3p . Using induction in the number of particles it is easy to show that the 
two contributions cancel each other, and therefore (|A.5j) is the final answer to the reflection 
factor. It also agrees with the results for the two-particle case in [2T] after accounting for the 
difference in the normalization conventions. 





A. 2 Reduction formula in the closed channel 

In the closed channel the Hilbert space is identical to that of the bulk theory, and the quantities 
of interest are the multi-particle matrix elements of the boundary state: 

{B\qi,qi; . . .;qN,qN)in = {2ti)^ 6C^qj)6(^qj)u} {qi,qi) {{qi,qi]) 

j j 

{uj{qi,qi) is a normalization factor introduced for convenience in order to conform with the 
conventions for and in the main text). In the bulk theory matrix elements are expressed 
in terms of the correlators via the reduction formula, but due to the presence of a boundary 
in time they need to be modified. First we express Af^{qi,qi) in terms of the in field as 

{B\Aliq,,q,)\q2, 92; . • • ; giv, qN)^n = / dydfe-'^^'i^''^^^^+''i^y-^''^^''dr (i?|a>,„(r, y, r)\in) 

V 2 j-00 

where I in) is a shorthand for \q2j q2] • • • j QNj QN)in- 

Now we use that for large negative r the 

interacting field approaches the asymptotic field 

lim ^>i„(r, y, r) = ^"^/^^(t, y, r) 

T—t — OO 

and also 
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to obtain the rule for the elimination of a particle from the initial state: 



{B\Af^{qi,qi)\q2,q2; ■ ■ ■■,qN,qN) 



l2Z 



oo /"O 

-oo J —oo 



[d^, - d^y - + - d{r){dr + iLo{qi,qi))} 



T,y,r)\m) 



Applying this procedure successively, the multi-particle matrix element can be expressed as 
the product of integro-differential operators 

r dy f dt j dfe-^^^^'^^^+^^s^+^^^ja^ -dl-[P + m^ - 6{T){dr + 

J —CO J ~oo J 



/2Z 



acting on the correlator. The reduction formula obtained this way differs from its bulk coun- 
terpart by the presence of the 6{t) term. It also differs from the analogous expression (|A.2p 
in the open channel by containing e~'^^^'^'^'^ instead of cos(A:x) and by the extra — i(5(r)a;(g, g) 
term. Despite these differences it can be shown that it also truncates the momentum space 
Green function and puts the momentum on-shell. To do so we rewrite the open channel one 
leg contribution to the Green function (IA.4P in terms of the closed channel 



r d^^^.rdk^_,,, dk 

J_^2tt 2tt (27r)«-i' 



-ikr 



iZ 



icr(m )dm 



^ — k'^ — k'^—m'^+it J2m ^2_^2_^2_^'2_|_jg 



In parallel with the open channel calculation we can eliminate the dependence on and k 
by performing explicitly the differentiations and integrations. Finally it remains to show that 
the operator 

dT e-^^('?'^> +a;(g,g)2 _ s{t) {iuj{q,q) + dr)} 



truncates the leg 

, 2tt 



iZ 



LO^ — uj{q, q)^ + ie 



+ 



ia{m)dm 



2m 



ijj^ — u{q , q)^ + ie 



q^ + vn? — m"^. We compute separately the contributions from + 



where, as before, q"^ 

uj{q,qf' and —5{T){iu>{q,q) +dr)- The second is simpler: it gives a factor —i{Lo{q,q) + u) 

and substitutes r = in the to integration. The operator + oj{q,q)^ gives the factor 
{—uj^ + uj{q, qf) and the r integration gives 



uj — uj{q, q) — ie 



-iV 



u}-ui(q,q) 



+ tt6{u! - uj{q, q)) 



Let us first concentrate on the a terms. The 5 function does not contribute due to the pre- 
factor —uj"^ +uj{q,qf'^ while the principal value can be replaced with the function itself and 
this just kills the a term coming from the 5{t) part. In the Z term we use that 

uj + uj{q,q) -i 



up' — u!{q, q)'^ + ie uj — u;{q, q) + ie 
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This kills the principal value term above and results in a factor two in front of the 6 (r) part, 
which together gives 

-iZB^^^y{-iq, iuj{q,q),q, . . . ) 

Combining with the pre-factor ^== in the reduction formula gives the contribution of one leg 
as 

In the closed channel we have just one type of disconnected parts, the one coming from the 
Green functions (|A.3p . which contains at least one bulk two-point function as a disconnected 
piece. Just as in the two-particle case this does not contribute via the reduction formula, and 
since it appears multiplicatively the whole contribution coming from the disconnected part is 
vanishing. As a consequence only the connected part of the Green function contributes to the 
boundary state, which is therefore identical to the q ^ iuj continuation of the result (IA.5I) for 
the reflection factor, and so K^~^^ can be expressed in terms of . We make this relation 
explicit for the two-particle term in subsection 3.3. 
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